Deformation parameter for diffuse density 
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In extracting deformation parameters from multipole moments for deformed nuclei, one commonly 
uses the formulas which are based on a sharp-cut density distribution. We discuss a possible am- 
biguity for this procedure and clarify the role of diffuseness parameter of the density distribution. 
For this purpose, we use a deformed Woods-Saxon density as well as a density distribution obtained 
from the self-consistent relativistic mean-field (RMF) model. We show that the formula using a 
root-mean-square radius instead of a sharp-cut radius requires a large correction due to a finite 
surface diffuseness parameter even for stable nuclei. An implication to neutron-rich nuclei is also 
discussed. 
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Deformation of a density distribution in the intrinsic 
frame is one of the most important concepts in nuclear 
physics. It is characterized by a finite value of intrinsic 
quadrupole moment Q, but a more intuitive quantity is a 
deformation parameter /3, which removes the trivial de- 
pendence on the radius of a nucleus from the quadrupole 
moment. In order to extract the deformation param- 
eter from the quadrupole moment, one often uses the 
formula which is obtained by assuming a sharp-cut den- 
sity distribution[ll 01 ■ To partly account for the devi- 
ation of density distribution from a sharp-cut function, 
sometimes one also uses a formula in which the sharp-cut 
radius is replaced by a root-mean-square radius (with a 
trivial constant factor) [1, 0, |a]. In general, the latter 
formula is believed to be better than the former when 
the surface diffuseness parameter of density distribution 
is large as in neutron-rich nuclei. 

In this paper, we systematically investigate whether 
that is the case. To this end, we first use a deformed 
Woods-Saxon density. This model has an advantage in 
that the deformation parameter is given as an input pa- 
rameter. This enables us to check whether the commonly 
used formulas lead to the correct value of deformation pa- 
rameter. We also derive the correction to the formulas 
due to a finite value of surface diffuseness parameter. We 
then discuss the deformation parameter of neutron-rich 
Mg isotopes using the self-consistent relativistic mean- 
field (RMF) model. 

Let us begin with the definition for the intrinsic 
quadrupole moment. 



Q = ^^Jdrpir)r'YM9). (1) 

When the density p{r) has a sharp edge, that is. 
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with po = 3A/{4:TtRq), A being the mass number of a 
nucleus, and 



the quadrupole moment Q is evaluated as. 
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to the first order of deformation parameter /3. One often 
takes i?o = 1.2 A^/'^ fm for the sharp-cut radius Rq Q]- 
For a sharp-cut density ^, the root-mean-square radius 
is calculated as 
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again to the leading order of (3. Therefore, the relation- 
ship between the quadrupole moment and the deforma- 
tion parameter given by Eq. Q can be also written as 
iil 
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The effect of finite surface diffuseness of density dis- 
tribution can be accounted for using a deformed Woods- 
Saxon density. 



p{r) = 
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In order to derive the correction term to Eqs. Q and 
©, we expand Eq. {T)) with respect to the deformation 
parameter (3 and keep only the first order term, that is. 



p(r)^po(r)-^i?o/3>^2o(^) 
dr 



where 
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R{0) = Ro{l+pY2m), 



(3) 



Notice that the angle dependence of the surface diffuse- 
ness parameter a |y, S 13 as well as the dependence of 
-Ro on the deformation parameter due to the volume con- 
servation 0] does not appear in Eq. (jHl since they are 



higher order terms of /3. With the density given by Eq. 
©, the quadrupole moment Q is calculated as, 
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Since dpo/dr has a finite value only in a small region near 
the nuclear surface r ~ i?o, we expand r^ around r = Rq. 
Following the same procedure as in Refs. [3, [l3 > we find 



Q-\I^R'oPop(l + 27r'^,, 



(11) 



to the order of (a/i?o)^- One can eliminate the depen- 
dence on pq in Eq. ()11|) using the condition for the nor- 
malization, 



A — dr p{r) 
This yields, 
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Furthermore, one can use the relation 
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to eliminate the dependence on i?Oj leading to 
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A similar formula can be found also in Ref. [g- 

We now investigate the performance of Eqs. Q, lO, 
and pSI) using realistic density distributions. To this end, 
we first compute the quadrupole moment |^ using a de- 
formed Woods-Saxon density, 0. Following Ref. |l5l |. 
we choose Ro = 1.31^^/^-0.84 fm in Eq. ©. This value 
was obtained by fitting to theoretical as well as to exper- 
imental density distributions for a number of nuclei with 
the Woods-Saxon shape '15'|. For a deformation parame- 
ter j3 in Eq. ©, we choose /3 — 0.3. Once the quadrupole 
moment is obtained, we can use Eqs. Q), ©, and H15() 
to obtain an approximate value for the deformation pa- 
rameter, 13. If the formulas worked perfectly, they would 
lead to /3 = 0.3 as in the original density distribution. 

Figure 1 shows the deformation parameters obtained 
in this way as a function of diffuseness parameter a in 
the density distribution. The upper panel is for A—AQ, 
while the lower panel for A=238. The dashed line is ob- 
tained with Eq. ||HJ) using the root-mean-square radius. 
We see that this formula significantly underestimates the 
deformation parameter, especially for the lighter system, 
A=40, except when the diffuseness parameter a is close 
to zero. Surprisingly, the formula does not seem to work 
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FIG. 1: Deformation parameters obtained with several for- 
mulas as a function of surface diffuseness parameter a in the 
density distribution. The deformation parameter is set to be 
13 = 0.3 in the deformed Woods-Saxon density. The upper 
panel is for A=40, while the lower panel for yl=238. The dot- 
dashed line is obtained with Eq. Q assuming a sharp-cut 
density, while the dashed line with Eq. using the root- 
mean-square radius. The solid line is obtained with Eq. (1151 , 
including the surface diffuseness correction. 



even for stable nuclei around a ~ 0.55 fm. When one 
takes into account the surface diffuseness correction with 
Eq. (|15|l . one obtains a reasonable agreement with the 
exact value of the deformation parameter, as shown by 
the solid line. The dot-dashed line is obtained by using 
Eq. igji with i?o = 1.2 A^/^ fm. We see that this formula 
provides a reasonable value of the deformation parameter 
if the surface diffuseness is around 0.55 fm as in stable 
nuclei, although the deviation from the exact value be- 
comes large when the surface diffuseness is around 1.0 
fm. 

Let us now discuss the deformation parameters ob- 
tained with self-consistent RMF calculations. The up- 
per panel of Fig. 2 shows the deformation parameters 
for Mg isotopes. The filled circles, the open circles, and 
the open triangles are obtained with Eqs. I0J, ©, and 
H15|l . respectively. In this calculation, we use the NL3 
parameter set |llj and assume the axial symmetry. We 
use the computer code RMFAXIAL 2] , which solves the 
RMF equations using the harmonic oscillator expansion 
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FIG. 2: Upper panel: Deformation parameters for Mg iso- 
topes extracted from the relativistic mean-field (RMF) den- 
sity distribution. The open triangles, the open circles, and 
the filled circles are obtained with Eqs. Eqs. Q, (|U, and 
H15|l . respectively. Lower panel: The diffuseness parameter 
for each nucleus obtained by fitting the RMF density to the 
Woods-Saxon shape. 



method. We employ the constant gap approach with the 
pairing gap given in Refs. |12lll3| . In order to estimate 
the diffuseness parameter a, we expand the density dis- 
tribution into multipoles (see Appendix B of Ref. [l4j). 

p{r,0) = Pair) + p2{r)P2{cose) + p4{r)P4{cose) + ■ ■ ■ , 

(16) 
and fit the monopole density po (r) with the Woods-Saxon 
shape, Q. The diffuseness parameter thus obtained is 
shown in the lower panel of Fig. 2. One finds that the 
deformation parameter estimated with Eq. jnj (the open 
circles) is considerably smaller than that estimated with 
Eq. ()15|l which includes the surface diffuseness effect (the 
open triangles), in accordance with the study with the 
deformed Woods-Saxon density. One also finds that the 
deformation parameters estimated with Eq. (@J) are close 
to those with Eq. H15|) (the filled circles), although it 
might simply be accidental. We have checked that our 
conclusions remain the same also for Sr isotopes. 

Figure 3 shows the potential energy surface for ''"Mg 
nucleus as a function of deformation parameter /?, ob- 



FIG. 3: Potential energy surface for Mg. The solid line, the 
dashed line, and the filled circles are obtained with Eqs. Q, 
@, and CSJi respectively. 



tained with constrained RMF calculations. The solid 
and the dashed lines are obtained with Eqs. Q and 
©, respectively. Since Eq. lO underestimates the de- 
formation parameter, the distance between the prolate 
and oblate minima is also underestimated. The energy 
surface which takes into account the surface diffuseness 
correction is denoted by the filled circles. As we showed 
in Fig. 2, it is similar to the energy surface obtained with 
Eq. ©. 

In attempting to extract the deformation parameter 
from an experimental quadrupole moment, an interest- 
ing question is how to estimate the surface diffuseness 
parameter a from experimental data. Unfortunately, the 
root-mean-square radius alone does not determine the 
radius and the surface diffuseness parameters, Rq and a, 
simultaneously. One possible prescription is to assume 
Rq = 1.31^1/3 - 0.84 fm as in Ref-O, and use Eq. 
(|14|l to estimate the surface diffuseness parameter a from 
the root-mean-square radius. In fact, our RMF calcu- 
lations show that the Rq parameter estimated from the 
fitting to the monopole density is well parametrized by 
this function up to the drip line nucleus, at least for the 
Mg isotopes presented in Fig. 2. One can then use Eq. 
(|13|l to estimate the deformation parameter /3 from the 
quadrupole moment Q. This procedure may be impor- 
tant in discussing the deformation parameter of neutron- 
rich nuclei, where the surface diffuseness parameter is 
expected to be significantly larger than that of stable 
nuclei [l6j. 

In summary, we discussed the role of surface diffuse- 
ness parameter of density distribution in converting the 
quadrupole moment to the deformation parameter. For 
this purpose, we used both the deformed Woods-Saxon 
and the RMF density distributions. We showed that the 
widely used linear order formula with root-mean-square 
radius significantly underestimates the deformation pa- 
rameter. After including the surface diffuseness correc- 



tion, the resultant deformation parameters were found to 
be close to those estimated with the linear order formula 
with a sharp-cut radius. 

The present consideration will be important when one 
discusses the deformation properties of neutron-rich nu- 
clei, where the surface diffuseness parameter is expected 
to be large. In particular, when one draws a two dimen- 
sional potential energy surface spanned by proton and 
neutron deformation parameters, /3p and /?„, it may look 
considerably different depending on which formula one 
employs in estimating the deformation parameters. Such 
studies are now in p rogr ess, and we will report on them 
in a separate paper |l7l |. Also, it will be an interesting 



future problem to extend the formula derived in this pa- 
per by including the higher order terms of deformation 
parameter /3. This will involve the angle dependent sur- 
face diffuseness parameter and the condition for volume 
conservation. 
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